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' Abstract 

^\ ' We study billiards in plane domains, with a perpendicular magnetic field and a potential. 

We give some results on periodic orbits, KAM tori and adiabatic invariants. We also prove 
the existence of bound states in a related scattering problem. 

> 

' Classical billiards are popular models for various physical systems, in fields ranging from me- 

chanics of systems with impacts and ergodic theory Q to semiclassical methods in quantum chaos. 
' In particular, billiards in a magnetic field appear to be relevant for the study of transport properties 

in mesoscopic systems, diamagnetism and the quantum Hall effect (see for instance Q). 
' In this work, we consider the classical motion of a charged particle in a plane domain, with a 

, perpendicular magnetic field of intensity B and a potential V{x.). We first discuss a method for 

proving existence of periodic and quasiperiodic orbits. Then we give some results on KAM tori 
and adiabatic invariants for billiards in a magnetic field. Finally, we consider the scattering on a 
' hard disc in crossed electromagnetic fields, where we prove the existence of bound states. Details 

Q , on the present results can be found in [|[ ^ § . 

. 

^ ■ 1 Billiards and Periodic Orbits 



O 

Bouncing Map: We consider the classical motion of a particle in a connected domain Q (which 
' may be unbounded or not simply connected). The boundary dQ is parametrized by its arclength, 

^ . x(s) ~ (A(s), Y{s)), with A'(s)^ + F'(s)^ — 1; the unit tangent vector and the curvature are given 

respectively by t(s) = (A'(s), r'(s)) and k(s) = X'{s)Y"{s) - X"{s)Y'is). 
Inside Q, the billiard flow is defined by the Lagrangian 

£(x, x) = imi^ ^ . _ y(^^^^ 

where A(x) = ^B{—y,x) is the vector potential in symmetric gauge (we will adopt the sign 
convention qB < 0). 

The dynamics is defined in the following way. Assume that the billiard particle starts on the 
boundary at x(so), with a velocity Xq making an angle with t(so)- It then evolves in Q according 
to the Lagrange equations. If the particle returns to the boundary, at a point x(si), with a velocity 
xi making an angle —9i with t(si), it is reflected elastically, meaning that it leaves the boundary 
again with an angle 9i (Fig. 1). The i-th collision may thus be parametrized by the coordinates 
{si,9i), or, alternatively, by the Birkhojf variables {si,Ui) = Zi, where u,; = x^ • t(si) = jxijcos^i 
denotes the tangent velocity. 

As long as the particle returns to the boundary, we may thus describe the dynamics by the 
bouncing map T : zq i— s- zi. 

*Procecdings of the 3^'^ International Summer School/Conference in Maribor, Slovenia, 24 June - 6 July, 1996. 
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Figure 1: The abscissa sq of the starting point and the angle 9q uniquely define the trajectory 7, 
and, if the boundary dQ is reached again, the coordinates (si, 9i) of the next collision. 



Generating function: A very useful tool to study the bouncing map is its generating function 
G{so,si) [p|. which satisfies 

dG dC 

= -uq, -—=ui. (2) 

OSo OSi 

It can be constructed in the following way: let 7 be a trajectory connecting the points x(so) 
and x(si), and F{so,si) — p • dx the action along 7, where p — d^C — mx + gA(x) is the 
momentum. We know from analytical mechanics that for infinitesimal variations of the end points 
dxi = t(si)dsi, the change of action is dF — — po • dxo + pi • dxi. It is then easy to check that 

G{so, si) = -F{so, si) + 1^ / ' Y{s)X'{s) - X{s)Y'{s) ds (3) 

satisfies the relations (|^). 

Here we have assumed that there is a unique trajectory connecting the two points on the 
boundary. In fact, there may be several or no solution of the Lagrange equations for given end 
points, so that the generating function may be multiply defined on some domain, and not exist 
on another one. This gives rise to new complications, but the main arguments presented in the 
following can be transposed to this more difhcult situation. 

Periodic orbits: The relation (||) is useful to compute periodic orbits. If sq, si, . . . Sn-i is a 
sequence of arclengths on the boundary, we define the n~point generating function 

G'(")(so, si, • ■ • Sn-i) = G(so, Si) + G{si,s2) + ■■■ + G(s„_i, So)- (4) 

The law of specular reflection implies that there is a periodic orbit of period n, hitting the boundary 
at x(so), . . . , x(s„_i), if and only if dG'^^^dsi — for each i (assuming G^"-' is defined and 
sufficiently differentiable) ; in other words, periodic orbits correspond to stationary points of G*-"-*. 
The linear stability of the orbit is determined by noting that 

dwo = -G2odso - Giidsi _ 5-^+^G . . 

dui - Gndso + Go2dsi '^^'^ ~ q,i^ Q^k ^ 
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implies dzi — r'(so, Si)dzo, where 

7"(so,si) = --^ ( „ „2 r ) 

Gil V G20(-f02 ~ Gil G02 / 

is the Jacobian matrix of the bouncing map. Since T' has unit determinant, T is area-preserving. 
For n iterations, dz„ = S'„dzo, where 5'n(so, . . . , s„_i) = T'(so, s„_i)T'(s„_i, s„_2) ■ ■ •T'(si, sq). 
The eigenvalues of Sn are A± = t ± vt^ — 1, where t = iTr^n. Hence, the orbit is hyperbolic if 
|ij > 1, elliptic if \t\ < 1 (under certain conditions, t may be related to second derivatives of G*^"^ 

. . . 

The center manifold theorem implies that hyperbolic orbits are unstable, even when nonlinear 
terms are taken into account. Elliptic orbits are generically nonlinear ly stable (in the sense of 
Lyapunov), as a consequence of the KAM theorem. A standard result is 
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Theorem 1 Let T be measure-preserving and in a neighborhood of the periodic orbit. Assume 
that the eigenvalues are such that (A^)"^ 7^ 1 and (A±)^ 7^ 1. There exists C, depending only on 
the second and third derivatives of T along the orbit (and which can be explicitly computed), such 
that if the non- degeneracy condition C ^ is satisfied, any point of the orbit has a neighborhood 
which is invariant under the map T". 

This result is proved by applying Moser's theorem to the Birkhoff normal form of T". 
2 Billiards in a Magnetic Field 

The particular case when there is only the magnetic field (F(x) = 0) was first considered by Robnik 
and Berry For a given energy E, the trajectories are arcs of Larmor radius /i = \j2raE j \qB\. 
The generating function and the Jacobian matrix can be explicitly expressed in terms of geometric 
properties of the boundary |^ . 

We first consider the billiard in a convex domain, i.e., we assume that the radius of curvature 
/9(s) = 1/k(s) is smooth and bounded by positive constants Pmin and Pmax- An important class of 
orbits are the "whispering gallery modes" : they correspond physically to quasiperiodic trajectories 
skipping along the boundary; in phase space, these solutions live on invariant curves which are 
close to = 0, TT. In the zero field case, existence of such orbits was proved by Lazutkin ||l^. As 
remarked in when a magnetic field is added, the dynamics near the boundary depend on the 
relative value of /i, pmin and Pmax- This is confirmed by the following result: 

Theorem 2 |^ Assume the boundary is and satisfies < Pmin < pis) < Pmax < 00. There 
exists a Cantor set of invariant curves with positive measure in the three following situations: 

1. For < /i < cxD, near ~ n. It corresponds to backward skipping trajectories, which are curved 
towards the boundary. 

2. For fj, > Pmax; near 9 = 0. It corresponds to forward skipping trajectories which are curved 
away from the boundary. 

3. For fi < pmin; near = 0. It corresponds to backward skipping trajectories starting with a 
forward glancing velocity, and performing almost complete Larmor circles. 

The proof relies on a perturbative expression of the bouncing map for small sin 6, which can 
by analyzed by Moser's theorem The difference between the regimes 0^0 and 6* ^ tt is due 
to the symmetry breaking effect of the magnetic field. When pmin < /i < Pmax, invariant curves 
near 9 = are absent due to discontinuities by tangency. Note that in contrast with a theorem by 
Mather p| , in a magnetic field invariant curves still exist when the curvature is allowed to vanish. 
In this respect, the magnetic field has a stabilyzing effect. 

This result can be extended to more general billiard domains. Consider for instance the bil- 
liard outside a given convex curve. One can try to construct a trajectory of the outside billiard 
by completing every arc of an inside trajectory to a full circle (Fig. 2). There is a one-to-one 
correspondence between inside and outside orbits if the following property is satisfied: 




Figure 2: Inside-outside duality. In (a), the /i-intersection property is satisfied, there is a one-to- 
one correspondence between inside and outside trajectories. If this property is not satisfied, the 
duality may be destroyed (b). 



3 



Definition 1 A closed plane curve is said to have the yu-intersection property for some > 
if any circle of radius fi intersects it at most twice. 



Lemma 1 ||3|] yl closed plane convex curve with extremal radii of curvature Pmin md Pmax satisfies 
the ^-intersection property if < Pmin or ^ > pmax- 

We conlude that the inside and outside biUiards are equivalent in low or high magnetic field. 
In fact, it is possible to show that the inside-outside duality remains true for backward skipping 
trajectories, even for intermediate magnetic field. Thus, Theorem ^ is valid for outside as well as 
for inside billiards, implying that this large class of billiards is not ergodic. 

More generally, one can consider domains which are not convex, but whose boundary has a 
bounded curvature: \k{s)\ < 1/pmin- It is possible to show that point 1. of Theorem ^remains 
true if p < Pmin) which is once again a manifestation of the stabilyzing effect of the magnetic field. 
In this case, it is of particular interest to consider the strong magnetic field limit. 

Proposition 1 |^ Assume the boundary is , fc > 3, and has a bounded curvature. For suffi- 
ciently small jjL, the bouncing map is C^~^ and takes the form 

,2 , 



si = So — 2/isin0o + M sin (?o a(so, 6*07 m) (mod 
9i = 6f) 



sin^ da b{sn, 00, fi). (7) 



The functions a E C'' ^ and b E C'^ ^ are uniformly bounded for s £ TR,0 < < n, \dQ\-periodic 
in s, and admit expansions in p which can be explicitly computed. 

The bouncing map (0) has the structure of a perturbed integrable map, where the factors sin^o 
ensure that the boundaries = 0, tt are fixed. One can again study invariant curves, of the form 
I{s,9) = const, using Moser's theorem, although the analysis is complicated by the fact that the 
frequency fl(9) = —2fism9 is multiplied by the small parameter p and is not monotonic. 

An alternative approach to KAM theory is to construct adiabatic invariants J(s, 9), such that 
J(si, 6*1) — J{so, 60) is as small as possible. 

Theorem 3 |^ If the boundary dQ isC'', k > 3, and has bounded curvature, there exists a function 
J{s,9) such that J{si,9i) — J{so,9o) +C'(p'^+^). If dQ is analytic, there exists a function J{s,9) 
such that J{si,9i) = J(so, 9o) + ©(e-^/'^l'^l). 



In fact, Theorem ^ is true for a large class of maps, including (^. In the case of a billiard, 

1 . . 2 



J(s, 9; 11) = 9 + fi sin ( 



■k(s) + ^p cos 6Ik(s)2 + e)(p2) 



(8) 



If the boundary is analytic. Theorem ^ implies that for any initial condition {sq,9q), {sn,9n) 
remains at a distance of order p from the level curve J(s, 9) = J{so, 9o), during a time of order e^. 



3 A Scattering System 

We now consider the case where the potential is given by a uniform electric field, ^(x) = —q£ ■ x, 
£ = (0, \£\). One can introduce dimcnsionless variables such that the Lagrangian becomes 

C = ^x^ + ^{yx-xy)-ey, (9) 

where e — \£m/qB^\ measures the strength of the electric field. 
The trajectories are cycloids of the form 

x{tp) = a + sip + pcos{^lj — ^p), y{^p) = b + psin{ip — 'ip), (10) 

where p reduces to the Larmor radius when e = 0. To account for the conservation of energy 
E ~ ^{e"^ + 2eb + p^), we introduce a second parameter p — \/2E — e^, so that p = a/p^ — 2e6. 
The range of the tangent velocity is |u| < p(l + 0(e)). 
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Figure 3: A trajectory scattered off the hard disc. Trajectories coming in from infinity leave the 
scatterer again with probabiHty one. However, some orbits may form "bound states" which are 
indefinitely bouncing on the scatterer. 

In m, the following problem was studied: consider the billiard outside a circular scatterer, 
parametrized by x(s) = (coss, sins). Are there trajectories which remain "trapped" in the vicinity 
of the scatterer, in spite of the drift due to the electric field? In fact, one can easily show that 
a particle drifting in from infinity will leave the scatterer again with probability one (Fig. 3). 
However, it is possible that some trajectories bounce on the scatterer indefinitely in the past and 
in the future, forming a classical "bound state". To show this, we first need to construct the 
generating function. Given two end points x(so) and x(si), we have to determine the equation of 
the corresponding cycloid (|lO|), and use the general formula (^). The perturbative result is: 

Proposition 2 |^ Let s± = (si ± so)/2. There are positive constants c\, ci and Eq, such that for 
ciE < s_ < TT — ciE, fi > 1 + ci£ and < e < Eq, the generating function of the bouncing map is 
unique, an analytic function of s±, fi and e, and admits the expansion 



where C , S , R, At/i denote functions of s_ alone: 

C = coss_, S' = sins_, R = [i^ - S^, ^ 2n - Atccos[1 - 2 / fi^]. (12) 

The point is that the generating function has the form G{so, si, s) = Go{si — So)+£Gi{so, si, e). 
If we substitute this expression in (^, we obtain that the bouncing map has the structure of a 
perturbed integrable map: 

Corollary 1 There is a positive C3 such that for \u\ < /i(l — c^e), fi > 1 + C2S and < e < Eq, the 
trajectory returns to the boundary, the bouncing map is analytic and of the form 



This result asserts that the particle will return to the scatterer if it starts with a sufficiently 
large normal velocity. The problem is now to show that some of these trajectories indefinitely 
return to the scatterer. This can be achieved once again by using Moser's theorem: indeed the 
existence of two distinct invariant curves in phase space implies the region between them to be 
invariant under the map. 

Theorem 4 g There is a positive Si such that for e < ei and /i > 1 + C2e, the scattering system 
has a set of bound states with positive measure. 

The problem with KAM theory is that one has in general very bad estimates on the bound 
£1. One can improve them by studying periodic orbits, which can be surrounded by islands of 
stability for much higher values of the electric field. As described in Section 1, we can compute the 
two-point function G*-^-* and analyze its stationnary points (for this purpose, we needed to know 
(O) at second order in e). In this way we can prove the existence of two orbits of period 2: 




AV^Msin^s+-ilAV^2 ^0(^3)^ 



(11) 



so + n{uo)+£f{so,uo,e) (mod 27r), 
uq + £g{so,uo,£). 



(13) 
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Figure 4: Orbits of period 2, for e = 0.1 and fj, = 1.5: (a) hyperbolic orbit, (b) elliptic orbit. The 
elliptic orbit is generically surrounded by a set of bound states. 



1. An orbit hitting the scatterer at s = 7r/2,37r/2, which is hyperbolic for small e (Fig. 4a). 

2. An orbit hitting the scatterer at s = (5, tt — (5, where 5 = sAip {it / 2) / R{tt / 2) + 0{e^), which is 
elliptic for < e < £2 (Fig- 4b). Using Theorem |l], we can show that if the orbit is elliptic, 
then it is stable with probability 1 with respect to d/ide. We have the estimation 

" 4 4AV'(7r/2) ^ ,^,2' 

which is in good agreement with numerics (roughly, bound states exist as long as the drift per 
cycle is smaller than the radius of the scatterer) . 

This method of searching elliptic orbits can be used to study existence of bound states for 
more general scatterers. In fact, if the billiard posseses elliptic orbits in zero electric field, they 
will generically survive small perturbations. Interesting open problems include (1) the existence 
of a critical electric field beyond which there are no bound states, and (2) an understanding of 
transport in phase space, and its influence on the transit time of a particle drifting in from infinity. 
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